Low temperature resonances in the electron 
heat capacity and its spectral distribution in 

finite systems. 

N.K. Kuzmenko^, V.M. Mikhajlov^ 

^V.G.Khlopin Radium Institute, 194021 St. -Petersburg, Russia 
^Institute of Physics St. -Petersburg State University 198904, Russia 

Abstract 

Temperature variations of the heat capacity (C) are studied in a low temperature 
regime T < dp ^ £f/N for 2D-, and 3-D-systems with N ~ 10^ lO'^ treated as a 
canonical ensemble of A^-noninteracting fermions. The analysis of C is performed by 
introducing function '^{e), the spectral distribution of C, that gives the contribution 
of each single-particle state to C . This function has two peaks divided by the energy 
interval Ae ~ (2 ^ 5)T. If at some temperature T^es there takes place a resonance 
i.e. the positions of these peaks coincide with energies of two levels nearest to ep 
then C vs T can show a local maximum at T^es- This gives possibility to assess the 
single-particle level spacings near the Fermi level. 
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1 Introduction 



The heat capacity of solids characterizes the temperature variation of the in- 
ternal energy corresponding to vibrations of lattice ions and movement of so 
called free electrons confined in solid. At temperatures T above the Debye 
temperature the phonon heat capacity prevails, fall in temperature below 
Td first equalizes the phonon and electron contributions to C and then at 
T ^ Tf) the main part of C belongs to electrons [T] . This general tendency in 
temperature evolution of C is observed in both macroscopic and mesoscopic 
systems though in the latter finite sizes can reveal themselves in the temper- 
ature and size dependence of T^j [2]. Inside the range, < T < Tjj, the 
relationship between the electron and phonon components is regulated by the 
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particle number N, the ratio of To to the electron Fermi energy (both quan- 
tities are material dependent). In mesoscopic systems the spatial shape and 
symmetry of a sample can also affect this relationship. However the phonon 
and electron contributions to C are practically independent that allows us to 
study here the low temperature electron heat capacity of mesoscopic systems 
separately from the phonon one. 

The electron movement in normal (nonsuperconducting) solids at low temper- 
atures can be described in the framework of the NONINTERACTING PAR- 
TICLE MODEL (NPM) provided particles are understood as quasiparticles 
introduced by the Fermi liquid theory f3]. This theory treats quasiparticles 
as compound fermion type excitations with an effective mass m* that differs 
from the free electron mass m and material dependence of m* causes the dif- 
ference in the Fermi energies ep of a solid and the Fermi gas in the same 
volume. These quantities, m* or ep, via po^ep), the electron state density in 
the vicinity ep, determine the prefactor in the linear in T Sommerfeld's law 
for the electron heat capacity 



in Eq. ([T]) and hereinafter T is measured in energy units. In macroscopic 
systems this law is fulfilled in a wide low temperature range, < T <^ ei?, ir- 
respective of the size and shape of the body. For some materials with m* ^ m 
there are experimental and theoretical evidences pointing out that the only 
parameter m* is insufficient to describe C in a wide temperature range since 
the heat capacity as a function of T can include a kink separating two linear 
temperature dependencies with different slopes [1]. The corresponding tem- 
perature is relatively high and we suppose that such radical change in m* 
does not concern a mesoscopic system at low temperature. In finite systems 
Eq. ([1]) is applicable down to an interval of T of order of interlevel spacings 
(~ Ep/N). Within it C increases exponentially and can display some local 
extrema which we shall consider in detail below. In mesoscopic systems there 
is also a transitional region of temperature quasilinearity where po{ep) de- 
pends on the shape of the body though these corrections are small ~ A^~^ and 



The temperature scale in finite fermion systems explicitly depends on ep that 
is material dependent: for electrons in metals ep ~ lO^-hlO^ K, in heterostruc- 
tures ~ 10^ K and for trapped Fermi gases ep ^ IfiK. Therefore, not to be 
bound to an individual mesoscopic system we measure all energies and the 
temperature in 6p, the average level spacing near ep 



C/ks = —pQ{ep)T, 




Ar-2 



6f = dp/poiep) ~ ep/N, 



(2) 
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dp is the average degeneracy near ep, for asymmetric systems as a rule dp = 2, 
for symmetric ones values dp > 2 will be specially pointed out below. Thus all 
results of our calculations in Sees. 3,4 can be projected on any real system if its 
^i? is known. Owing to Eq. ([2]) the term "low temperatures" acquires a more 
definite meaning. In the first instance it will be used for the temperature range 
smaller than or of order of 6p but these temperatures have to be really low, 
much less than To otherwise the electron heat capacity cannot be observed. 

In the finite Fermi systems [7j besides the effective mass (we suppose that 
m* is of the same order as in macroscopic bodies) the quasiparticle descrip- 
tion makes use of several additional parameters pertaining to the quasiparticle 
mean field (the depth, radius, diffusivity of the potential and others). These 
parameters are formed mainly by the ion surroundings and among real meso- 
scopic systems there are cases of potentials close to rectangular 3D-wells as 
it occurs in metal clusters p] whereas in quantum dots the best description 
of electron shells is given by 2D-oscillators [9]. These cases will be used in 
illustrative examples of Sees. 3,4. Such description idealizes and simplifies the 
real picture of the electron movement in mesoscopic systems as it neglects ef- 
fects of disorder. Coulomb and phonon interactions. However we believe that 
detailed ascertainment of the mean field role can help to reveal the impact of 
the other factors on measured thermodynamic quantities. 

Independently of the type of electron confinement and of that whether a sys- 
tem is integrable or not the common feature of quantum energy spectra in 
different finite systems consists in the appearance of such irregularities in the 
eigenvalue distributions as level bunchings or high degenerated levels in sys- 
tem possessing some symmetry. These irregularities can reveal themselves in 
measured quantities e.g. in the low temperature electron heat capacity. This 
fact was for the first time established by Frolich [lOj for a model spectrum with 
high degenerated levels. Later an analogous result was obtained for spherical 
alkali clusters by Brack et al [11] where the single-particle potential was found 
by the temperature Hartree-Fock method. In both cases the low temperature 
variation of the heat capacity displays local maxima at those particle num- 
bers which correspond to completely occupied high degenerated levels. On the 
other hand the calculations of Denton et al [12] showed that a uniform single- 
particle spectrum with spin degenerated levels and equal spacings between 
them practically gives a monotonous increase of C if the particle number con- 
servation is exactly taken into account. The comparison of these results has 
stimulated us to look for the cause of such different behavior of C vs T and 
the investigation of this problem is the subject of the following sections. 

Our approach to this problem supposes that each mesoscopic system possesses 
a quite definite shape, linear size and particle number. Now such systems can 
be fabricated by using the modern technologies and examples of such possi- 
bilities are described in Refs. [nill3fl4fl5lll6fl7] . Thus our description of meso- 



3 




scopic systems follows the line of Refs. |10|[TT] but differs from the approach 

developed by Kubo p^l9ll20|[2T] . This approach takes into account that in 
many cases shapes, sizes and particle numbers of mesoscopic systems cannot 
be fixed exactly and so thermodynamic quantities of such systems have to 
be calculated with single -particle spectra defined by means of statistical dis- 
tributions. However even in that case when experimental conditions make it 
possible to study only a set of systems parameters of which are given in some 
limits our results can be employed through an averaging procedure. 

The particle number conservation plays an important role in low temperature 
investigations of a mesoscopic system. As shown in Refs. [TTIT2] the results 
of the canonical and grand canonical calculations can noticeably differ in this 
temperature range. In Ref. [2^ we have developed a canonical polynomial 
method appropriate for studying finite systems in wide ranges of T and A^. 
Here the method is modified to give a convenient qualitative description of 
the heat capacity and show the resonance character of amplifying the heat 
capacity at some temperatures in the low temperature range. 

The subjects mentioned above are distributed over the following sections. In 
Sec. 2 we introduce a new representation for the canonical occupation numbers 
as a continuous function of energy and then on this basis we construct the 
spectral heat capacity (p{e). The latter is especially appropriate to analyze the 
influence of the irregularities in the single-particle spectrum on C in NPM 
that allows us to explain in Sec. 3 the cause of arising the low temperature 
maxima in C. In this section several examples of such maxima are given for 
both asymmetric and symmetric systems. In Sec. 4 we consider variations of C 
caused by the shape variations in the 2D-, and 3D-mesoscopic systems at an 
invariable area or volume i.e. at the particle number conservation. Conclusions 
are given in Sec. 5. 



2 The spectral heat capacity distribution 

Thermodynamic properties of the CANONICAL ENSEMBLE (CE) of A^ in- 
dependent fermions are traditionally treated in the framework of the approx- 
imate method which virtually proceeds from the assumption that CE occu- 
pation numbers (n^) of single-particle states (their energies and degeneracies 
are Eg and ds respectively) are the same as in the GRAND CANONICAL 
ENSEMBLE {GCE) i.e. Ug = fs where fs is the Fermi-Dirac function 



1 + e'^^^-") 



P = T 



-1 



(3) 
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In Eq. ([3]) A is the chemical potential chosen so that the ensemble averaged 
particle number at each T is equal to a fixed value of in a studied system 

oo 

N = J p{e)f{e)de. (4) 





In Eq. (jlj) and hereafter we suppose that single-particle energies e > and 
p{e) is the level density 

p{e) = ^5{e-e,)ds, (5) 



which depends on N via Sg due to the connection between the size and particle 
number in mesoscopic systems. 

In the exact canonical description the occupation numbers Ug are such that 
the equation analogous with Eq. (jl]) is provided automatically 

oo 

N = j p{e)n{e)de. (6) 





irrespective of the temperature and of the energy count point whose role in 
GCE is played by A. 

T and A are treated as independent variables in the classic definition of GCE 
while in CE those are T and N (in both cases the volume ~ keeps fixed 
therefore below the values of the heat capacity C will be considered only at 
this condition, C = Cy)- In order to stress the temperature dependence of A 
at the description with rig = fg this approximate method will be designated 
by the term the EQUIVALENT GCE {EGCE) [23]. 

In the EGCE formalism the partition function Zegce can be constructed so as 
to describe occupation numbers with help of fs and simultaneously reproduce 
the correct relation between the free energy Tegce and the thermodynamical 
potential f2 

n = -T\TyZGGE\ ZGCE = ViS^^exv\-ii{eg-\)\Y', (7) 

s 

Tegce = AA^ + 1] = -T In Zegce (8) 
Eqs. d?]) and ([HD can be consistent if 

Zegce = exp{-j3\N)ZGCE (9) 
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that immediately leads to {ns)EGCE — fs- 

Therefore in the EGCE description the internal energy Eegce and heat ca- 
pacity Cegce are represented in the standard form 



d f 
Eegce = -^In^EGCfi = j p{e)ef{e)de, (10) 



dp 



OEegce _ ^2 



dT 



dp' 



IuZegce ^ j P{^)^ 



df{e) 



dp 



de.{ll) 



The temperature variations of A in the EGCE theory explicitly manifest them- 
selves in df /dp which depends on pdX/dp 



dm_(^ , ,dX\dfie) 
dp [ ^dpj Pde ■ 



(12) 



The quantity pdX/dp is straightforwardly calculated allowing for the temper- 
ature independence of N in the EGCE theory 



dN 
'dp 



dp 



pde 



pde 



(13) 



The first integral in this equation is the particle number dispersion in both 
GCE and EGCE 



dm 

pde 



de = {N') - ({N))' 



(14) 



{N^) and are the ensemble averages of operators iV^ and N. The 

quantity (A -|- pdX/dP) defines also values of dEEGCE/dN, if N and T are 
high enough and dE/dN makes sense (at small variations of N the single- 
particle spectrum is fixed for an N particle system on account of the volume 
conservation) 



d Eegce _ . . n^X 
~d^-^^^dp 



(15) 



Hence one can obtain that {S is the entropy) 

dJ-'EGCE 



dN 



dS _ dX 
' dN ~ ~df 



(16) 
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In consequence of Eqs. ( !T2|) and (fT3i) the integral for Cegce can be represented 
in such form that the integrand is obviously positive 



CEGCE = f3' J p{e)el6-X-P 

^ 



am 



de 



/3de 

f{e){l-f{e))de. (17) 



Below to emphasize the role of the term pdX/dp, we adduce results obtained 
without this term {f3dX/df3 = 0). Such results will be denoted by letters GCE. 

The exact canonical partition function in the polynomial representation was 
suggested in Ref. [22] • Here, to have a uniform description for EGCE and CE 
we introduce a factor exp {—(3XcN) into the partition function of Ref. [22] : 



Z{N) = [N] exp i-pXcN) 



where designations [A^] and [n], that appears below, are employed for sym- 
metric polynomials the orders of which are N and n respectively. These poly- 
nomials are defined in a space of variables g<j 

qs = exp{-f3{e, - A^)}, (19) 



e.g. 

[1] = ^^qs] [2] = ds^qs^ds^Qs^ + ^^(^« ~ 

S Sl<S2 s 

Recurrent relations for polynomials and convergence of series for [A^] are con- 
sidered in Ref. [221. 



The auxiliary parameter Ac in Eqs. (1180 . (fT9i) separates the energy regions 
above and below the Fermi energy Sp and serves as a point to count single- 
particle level energies. However, Z{N) in Eq. (fT8|) does not comprise Ac at all 
and this independence from Ac is inherent in all functions generated by Z{N), 
i.e. 

Z{N) = [N] U.=o . (20) 

As shown in Ref. [22] the polynomial [A^] can be expanded in a rapidly con- 
vergent series in products of n order polynomials, n < A^. In these series each 
term corresponds to n-particle-n-hole excitations with respect to the ground 
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state at T = in which all single-particle states are filled up to and including 
the Fermi level (F). The simplest example of such expansion is Z{N) for a 
system with np = dp {np is the occupation number of the Fermi level at 
T = 0). 

Z{N) = Zo{N)Z{N), 
Zo{N) = exp {-f3Eo{N)) ; Eo{N) = J] 

s<F 

N 

Z{N) = 1 + Y.[n][n], (21) 

n=l 

[n], [n] are symmetric n-th order polynomials. The particle polynomial [n] is 
composed of (e^ > K), Eq- (121]) . the hole polynomial of {Sg < Ac), 

Bp < Xc< ep+i, (22) 



F + 1 or in the general case F + k and F — k are the first or A;-th level above 
and under F respectively. If F is incompletely occupied at T = or there 
is a necessity to count energies not from Ac given by Eq. (!22|) but from ep+i 
or £F-fc then single-particle energies have to be summed in Eq. (I?T]) up to 
{F + i)-th or only down to {F — k)-th level whereas the series in particle 
and hole polynomials have to begin now with [0][fi] and [z^][0] respectively, 
[0] = [0] = 1; /i = J2s<F+i ds — N; u = N ~ J^F-k ds- By using this possibility 
one can write down the partition functions of N ± u particles as following 



Z{N ±u)= exp{-/5 [Eo{N) ± Ac(iV)z/]} { ^"''^ ^ (23) 

where Ac(A^) is chosen so that Ac > Bp{N) i.e. near ep for the A^-particle 
system and Eo{N) = J2s<f{n) ds£s, -^(^) is the Fermi level of particles. 

To find the canonical occupation number 

Us = -^^lnZ{N) = ^^lnZ{N) (24) 



(ug characterizes the occupation probability of level s for only one electron, 
the full number of particles in state s is ngds) 

oo 

N = J p{e)n{e)de = ^ dgUs (25) 
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we take advantage of properties of symmetric polynomials 



a d ^'^-^ 

J^\N] = - E i-T"-' \N + m] = (26) 

"■s ys m=l 

N 

= E (-)"'^' - ^] C (27) 

m=l 

In Eq. (1261) M is the maximum polynomial order for the single-particle basis 
we take into account, as a rule M ^ N . Both Eqs. (126|) and (!27|) are exact. 
Thus any of them can be used to determine n^. Nevertheless the convergence 
of the series can be hastened if Eq. (l26l) is applied for states with e < Sp 
(g^^ < 1) and Eq. (j^ZD for states with e > Sp {qs <l) 



M-N 

n{e < £f) = 1 - E (-)"+' exp[/3m(£ - A„+)], (28) 

m=l 

N 

n{e > £f) = E (-)'"^' exp[-/3m(£ - A„_)], (29) 

m=l 

exp[T/?mA„±] = Z{N ± m)/Z(iV), (30) 
m\m± = ± [^(iV ± m) - ^(A^)] . (31) 

Eqs. (130|) and (13T!) indicate that each power m in Eqs. (128|) . (129!) enters into 
?2(£:) with its own chemical potential not arbitrary chosen but conditioned by 
the difference of the canonical free energies, Eq. (I3T!) . which have to be calcu- 
lated at invariable volume, i.e. with the same single-particle energy spectrum 
as in the case of particles. 

The comparison of Eqs. fl28l) and fl29l) with similar expansion of f{e) 



f{e < ^f) = 1 - E i-r^' exp[/5m(e - A)] (32) 

m=l 

oo 

f{e > £f) = E ^-T^' exp[-/3m(£ - A)] (33) 

m=l 

shows that the distinction between n{e) and f{e) is caused in general by the 
difference between Am± in Eqs. fl28l) . fl29l) and the GCii^ chemical potential 
A, though at high enough N and T ( the case considered in Eq. ([T6|l ) they 
become indistinguishable at least for small m <^ N when 
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The difference between A and Xm± can be assessed through the relation of 
Zegce and Zqe obtained by expanding Zqce, Eq- ©, in symmetrical poly- 
nomials 

N M~N 

Zgce = [iV](,) +EW-m\x) + E W + m\x) ■ (34) 

m=l m=l 



The subscript (A) is applied in Eq. to point out that polynomials [A^ ± m]^^^ 
are composed of Qs, Eq. (fT9|) . in which Ac is replaced by the GCE potential 
A. Hence, Zqce gains the form including potentials Xm± defined by Eqs. (130|1 . 

m- 



Zegce = e-^^"" Zgce = Z{N)R{N); 

N M-N 



m=l m=l 



As shown in Eqs. ffTS]) . fl2Ul) . Z(A^ ± m) = Zge{N ± m) does not comprise A 
or Ac. Thus, independently of values A and \m± the quantity In R{N) > and 

^EGCE < ^CE (36) 



At T — > the EGCE and Ci? free energies are identical. Consequently 
Eq. (135|) for superlow temperatures gives an estimation: 

Am+ > A > Am-- (37) 



In another temperature regime, 5 <^ T <^ sf (5 is the mean level spacing), 
the ratio Zegce/ Zce ~ (^/^)^^^; Refs. [T2j, [22], i-e. increases with T. It 
is possible if potentials A, Xm± comply with Eq. (157|) . At last, if a system 
is heated up to T > ^i? and survives at such temperatures both description 
(canonical and equivalent grand canonical) lead to the Maxwell distribution, 
i.e. n{e) = f{e), that can be provided if Zegce{N)/Zce{N) is a temperature 
independent constant. In this case, at T > Ef and A^ m it can be easy 
found that 

T 

Am+ ^ Am_ ^ A = -7 In — , (38) 

Ef 
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7 characterizes the increase of the average level density with energy 
/ e \ '^~^ 

Pa.(£) = 7— — , (39) 

e.g. 7 = 3; 2; 1 for 3D-, 2D- and ID-oscillator and 7 = 3/2; 1; 1/2 for 3D-, 2D- 
and ID-rectangular potentials. Thus, the values of the chemical potentials, 
Eq. fl38l) . are again in accordance with Eq. fl37j) . Therefore we can infer that 
Eq. fl37|) is justified in a wide temperature range. 

The estimations given by Eq. fl37j) lead to the well known fact( [24J) that at 
low temperatures n{e < ep) > f{s < ep) and n{e > ep) < f{e > Sp) i.e. 
the stepwise variation of the canonical occupation numbers {n{e)) vs e near 
Bp is pronounced more distinctly as compared with the grand canonical ones 

The temperature derivative of n{e) i.e. dn{e)/dT = —P'^dn{e)/dp is formally 
similar to df{e)/dT if the latter is expanded in powers of exp [±/3(£ — A)], 
Eqs. fl32|) and fl33|) . The distinction consists again in the chemical potentials: 
Xm± in n{e) and the only value of A in f{e), Eq. f|T2l) . 



AI-N 



^2 j2 (_)™+i^ U-A„+-/5 



dXr 



m=l 



dp 



dn{e < Ep) 
df 

exp [Pm{e - A.m_, 



(40) 



TV 



\m+l 



m \ 6 — Xr 



dXm- 

d/3 



dn{e > Ep) 
df " 

exp [-Pm{e - Xm+)] 



(41) 



m (^Am± + (3^^ = T-^ [In Z{N ± m) - In Z{N)] = 

= ±[E{N ±Tn) - E{N)]. (42) 

Eq. (H2|) indicates that for high N and T (when dE/dN makes sense) 



'"^^^ d(3 ^ ON 
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that is in accordance with Eqs. (fT5i) and (fT6|) . 



The relationship between the chemical potentials Xm,± and A establishes the 
general tendency for differences Eegce — Ece and Segce — Sce since function 
R in Eq. (1531) owing to Eq. (I57|) is an increasing function of T (_R — > 1 at 
T — > 0; -R — ^ M + 1 at T — oo) or a decreasing function of /5 = . This 
gives rise to 



E 



EGCE 



-—\^Zegce-Ece-^-^>Ece 



OR ^ 



g/3m{A-A™_) _|_ 



m=l 



(A„.+%)-(A + 4) 



-/3m(Am+-Am) ^ g 



Segce — Sce + Ini? — ^ttt > 'S'ce 

it ap 



(43) 



(44) 



As at low T only several first terms in Eq. (jUl) are essential the positive 
definiteness of {—dR/d(3) results in 



A„. + /3-^>A + ;)^>A„_ + /3-^ 



(45) 



at least for small values of m and low T. 

For calculations of the heat capacity C and for qualitative explanation of 
temperature and particle number dependencies of C at low T it is expedient 
to introduce into consideration the spectral heat capacity distribution ip{e) 
which at 6 = Eg gives the contribution of one electron on level s to C 



C = Y.ds'^i^s) = / p{e)ip{e)de 

s •' 

dn{e) 



^{e) = (e-eo)- 



dT 



(46) 
(47) 



An arbitrary constant Eq does not affect the values of C as 



/p(.)^& = 

(48) 

For convenience' sake of comparison with EGCE calculations below we keep 
Eq to be equal to (X+pdX/dP). An analogous spectral heat capacity for EGCE 
can be introduced on the basis of Eq. (ITTI) 
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The main feature of both function (p{e) and ipeff{s) important to interpreting 
C at low T consists in that these functions being depicted vs energy turn 
out to be two humped curves. Such function dependence on e can be easily 
exposed for ipeff{e) that is given analytically, Eq. (H9!) . and determined by 
only two parameters A and dX/df3 (the latter at low T is much smaller than 
A). In fact such form form of the curve for ipeff{e) arises due to the crossings 
of the bell-wise curve f{e) [1 — f{s)] with the parabola {e — X — l3dX/ d^Y and 
two maxima of {peff{e) at [3dX/d[3 <^ X are in points 



E^:,ii\T)^X±2AT (50) 
i.e. the energy distance between the maxima is 

AE(:£)(T)^4.8T. (51) 

Heating moves apart these maxima and thereby extends the energy space 
under the curve of ^{s). 

The explicit form of ^{s) depends on the single-electron spectrum structure 
more essentially through the set of parameters Xm± and dXm±/dj3. Neverthe- 
less ip{e) is similar to ipeffi^) on the whole though the curve for ip{e) passes 
under the curve for ipeff{e) at a fixed T on account of the relations between 
A and Xm±, Eq- (!37|) . The distinctions between ip{e) and ipeff{e) are discussed 
in the following sections. 



3 Origin of the low temperature resonance in C 



The two humped shape of ip{e) and ipeff{e) suggests an idea how the ap- 
pearance of the low temperature maximum in the temperature dependence 
of C (mentioned in Introduction) could be explained. At rather low temper- 
atures (T < Sp) both peaks of the spectral heat capacities can be so sharp 
that just under each of them there can be only one single-electron level (reso- 
nance levels). Sometimes a group of levels close together can take the role of 
the resonance level. In systems with even two such levels could be levels F 
(completely filled at T = 0) and F+1. In this case X+/3dX/d/3 ~ {eF+£F+i)/2 
and a local maximum in the temperature dependence of C could appear if the 
positions of two maxima of ip{e) or ipeffi^) coincide with the energies of these 
levels i.e. F and F + 1. The case when the maximum in C can arise at the 
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coincidence of only one peak of <f{e) with one level is given by such spec- 
trum near ep in which the level F is half-filled at T = and two levels F 
and F + 1 are remote from other single-particle levels of the spectrum. Then 
A -|- jSdX/dp — Ep and an amplification of C will be caused by the coincidence 
of the right peak of ip{e) with the position of level F + 1. 

The distance between the maxima in the spectral heat capacity depends on 
the temperature. For {peff{e), as shown in the previous section, Eq. flSTl) . it is ~ 
4.8T, for the canonical function ip{e) this distance is smaller, (2T-h4.8T), that 
will be shown in this section. Therefore only at a quite definite temperature 
that can be called the resonance temperature (T^es) the positions either of two 
peaks of ip{e) or one of them and the minimum of ip{e) exactly correspond to 
energies of the electron levels. Hence it follows that Tres is of order of some 
part (0.2 -h 0.5) of Ae, the spacing between the resonance levels. 

The condition T = T^es is necessary but insufficient for the appearance of the 
local maximum in C as it does not provide yet small values of the spectral 
heat capacity {(p{e) or ipeff{e)) at energies of single particle levels adjacent 
to the resonance ones. Indeed, if the peaks of (p{e) are not too narrow and 
other levels (with energies e^) divorced from resonance ones by spacings ~ 
Aep+i^p = Ep^i — Sp give noticeable values of ^{Ss) as compared with values 
of '{'{sp) and 9?(£f+i) then rise in temperature above T^es does not decrease the 
heat capacity and a local maximum in C will be practical invisible. Though 
both the positions of peaks and their widths are determined by T (increasing 
the temperature causes moving apart and broadening the peaks) the possibility 
of the local maximum appearance depends not only on T but also on the level 
structure near sp. 

Mesoscopic single-particle spectra are in general heterogeneous and local con- 
centrations of levels (level bunchings) alternate with their rarefactions. How- 
ever in some cases a spectrum near sp can be considered practically uniform 
with equal spacings between levels (the ideal equal level spacing spectrum with 
level degeneracies d = 2 is the spectrum of the ID-oscillator). For such spec- 
trum with d = 2 the canonical partition functions are given in Refs. [T^|E2] . 
In the case of the completely occupied Fermi level at T = these functions 
lead to the following equations for Xm±, if only terms greater than are taken 
into account (g = exp{—[35), 5 is the level spacing): 

Ai+ = ^i? -|- 5 — Tina, Ai_ = ^i? -|- Tin a, 

(52) 

A2+ = £p -\- 5, A2- = Sp, 

a = 2(l + g2)(l + 2g)-\ 
X = ep + 5/2 = {ep + ep+^)/2. (53) 
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It is easy to check that Xm±, Xm± + Pd\m±/dp, (m = 1,2) obtained in 
Eqs. (I52D, ([53D for T < 5 are satisfied Eqs. §7^), (gS]) even at T = 5. 

Since Ai+ ^ A2+ ( the difference is small at low temperatures) and terms with 
Am+ {m > 3) are inessential at T < 5 function (p{e) in the region e < ep can 
be approximated by ipeff{s), Eq. poj) . with A = Ai+. It means that the left 
peak of ip{e < ep) will be ~ 2.4T distant from the point ep+i — Ai+ ~ + 5 
(2.4T is the distance between the minimum and maximum of ipeff{s) that gives 
4.8T in Eq. flHT]) ). Thus, the maximum of ip{e < Sp) is closer to the point A, 
Eq. fl55]) . than the peak of ipeff{e) (the latter is 2.4T distant from A < Ai_|_). 
The analogous arguments point out that the right peak of (f{e > sp) is also 
closer to the point A than the right peak of (peff{s) at the same temperature. 
Thus, the peaks of ip{e) are roughly two times closer each other than the 
peaks of ipeff{e). Hence the positions of the peaks of ip{e) will coincide with 
Ep and e_F+i = £p + 5 at T^es — 0.35 > as it displayed in Figs. [H [2J 

Fig. [2] shows in addition that under the curve of ^{e) there are two more levels 
with noticeable (as compared with (^{ep) = ip{ep+i)) values of (p{ep-i) and 
ip{ep^2)- As mentioned above, it implies that the maximum in the canonical 
heat capacity will be practically absent while in the EGCE heat capacity 
it is well seen. Fig. [TJ For odd systems with uniform spectra near ep the 
temperature Tres turns out to be even more than that for even ones as the 
chemical potential A in this case ~ ej? while Xm± are practically the same as 
in Eqs. (!52l) (!53|) . i.e. the role of the resonance levels is performed by levels 
F — 1 and F + 1 that increases Tres, broadens the peaks of ip{e) and so the 
maxima in C{CE) and C{EGCE) are smoothed out. 

In Fig. [Hand some following figures we give instead of C vs T values of C /Cun 
to more distinctly display appearing the local maximum in C since in some 
cases such maximum looks more like a bend than a real maximum as it take 
place for C{EGCE) vs T in Fig. [T^. Cun and pQ{ep) are given by Eqs. ([T]) and 
fl5^ respectively. 

Looking at Fig. [2] one can assume that increasing a gap between levels F and 
F — 1 and simultaneously between F + 1 and F + 2 reduces contributions of 
levels F — 1 and F + 2 to C and so the maximum in C can arise. This is 
confirmed by calculations of C for a system with two isolated levels {F and 
F + 1) where F is completely occupied at T = (with np = dp = 2), Fig. [31 
In this case Tres — 0.3(£:i?+i — ep). For the same system with one electron on 
F ai T = {np = 1) the maximum also exists but Tres — ^-^^isp+i — £p) 
because now only the right peak of ip{e) coincides with level F + 1 and the 
minimum of ip{e) practically falls on the point ep, Fig. [H 

The example of a three level bunching generating local maxima in C for both 
even and odd systems is given in Fig. [5] The single-electron spectrum produc- 
ing the heat capacity variations and given in this figure belongs to a briquette- 
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shaped cavity with hard walls and with such ratio of the lateral lengths that 
removes accidental degeneracies i.e. each level is only spin degenerated {d = 2). 



Thus, in asymmetric system having as a rule only spin degeneracies of levels 
the low temperature maximum in C serve as evidence of a level concentration 
near F and the resonance temperatures (at which the maxima is observed) 
give information concerning level spacing Ae {Ae can be equal to e_F+i — 

or Ep+i — Ep^i) 

Tres ~ (0.2 ^ 0.4)Ae. (54) 



Additional investigation of particle number variations of C can help to ascer- 
tain detailed structure of the single-particle spectra. 

In the single-particle spectra of symmetric 21)-, and 3D-systems the majority 
of states is high degenerated. Such extreme concentrations of levels affects 
the low temperature variation of C via two factors both of which promote 
the appearance of the local maximum in C. The first is the extension, in 
the average, of spacings between high degenerated levels in comparison with 
asymmetric systems. The second is the approach of the canonical chemical 
potentials \m±^ Eqs. fl30l). fl3Tl) . to the grand canonical value A with increasing 
degeneracy d (i.e. increasing A^) that decreases the difference between (p{e) 
and ipeffi^) and makes more probable the appearance of the maxima in C as 
the peaks in (peffi^) are rather narrow. At very high d the canonical potentials 
become equal to A that results in identifying the canonical and effective grand 
canonical values of !f{e) and C. 

As an illustration of these properties of symmetric systems we will consider a 
model example in which two solitary equally degenerated levels F and F+1 are 
divided by the spacing 6 and at T = the level F is engaged (n^ = dp ^ 1). 

The polynomial method [22j allows the ratio of the partition functions, Eqs. fl30l) . 
to be written through finite sums in powers of q: 



Z{N±m) 
Z{N) 



Rm exp <^ =F/? 



q = exp(-/3(5); 



(ef + -) ± - 
^ 2^ 2 



(55) 
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Hence, the canonical chemical potentials Xm± gain the form 

Xm± = A ± 



6 6 
l-2{f35m)-HnR^]-, X = eF + -. (56) 



Since at d 1 values of Rm are proportional to d"^ the potentials Xm± ex- 
plicitly depend on Ind i.e. as it is followed from Eq. fl56|) the more d the closer 
Ami to A. However the general tendency given by Eq. fl37|) is not broken at all 
variations of d, that we demonstrate for two cases: d^q < 1 and d'^q ^ 1. 

If degeneracy d is not too high 

K ci < exp(/55/2); \nd< 136/2 (57) 

then neglecting terms ~ d^^ and introducing a new parameter p 

p = d^exp{-p5/2) <1 (58) 

one can retain in the sums for Rm only terms < that simplifies series for it: 



B,„ = d"'(mlr„)-' (59) 



Thereby Eq. fl37j) is satisfied as 

A„± = A±(a + 6)-; (61) 

a = 1-2{PS)-Hnd>0; 
b = 2{(36k)-^[ln{k\) +lnrk], 

the positive definiteness of a and b is the consequence of Eqs. (IFHl and (1^ . 

In the case of very high degeneracies 

d2exp(-/55) > 1 (62) 

series for Zm in Eq. (l55l) can be reduced (again neglecting terms ~ d^^ <^ 1) 
to the Bessel function of purely imaginary argument Im{y), y = 2dq^^'^: 

Zm ^ q-'^'^Uy) = E mm + k)\r' f (63) 
k=o ^^-^ 
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Now recollecting that the Bessel function asymptotic behavior {y ^ 1) does 
not depend on m one obtains by using Eq. (!55l) 

Rm = exp{(3Sm/2) (64) 



that implies on account of Eq. flSBl) the identity of Xm± and A. Thus, in the 
case of very high degeneracies (at >> 1) the canonical and effective grand 
canonical descriptions give the identical values of (f{e) and C. 

Here it is worth mentioning one more detail distinguishing symmetric sys- 
tems with d ^ 1 from asymmetric ones with d = 2. In the latter Cce near 
Tres is always less than Cegce due to the relationship between the chemical 
potentials \m± and A, Eq. ( 1371) . However in systems with high level degen- 
eracy approaching \m± to A can cause a so small difference between these 
potentials that the pre-exponent factors in C can compensate this difference 
and make Cce > Cegce in a narrow temperature range around Tjg""). E.g. 
in low temperature approximation the just considered two level model with 
d'^ exp(— < 1 and d = np ^ 1 gives: 

Cce ^ m'd'exp{-p6y, Cegce ^ m\d/2) exp{-p6). (65) 



For d = Hp = 10, e.g., it means that Cce > CpccE in the vicinity of 
Tres — 0.2 6 as it is shown in Fig. [61 

The local maximum appearance in C in symmetric systems strongly depends 
on that to what extent the Fermi level is occupied at T = 0, Fig. [61 Since in 
such systems the spectral heat capacities (p{e) and ipeff{e) are rather close, 
therefore, to simplify the problem we will address the later as it is characterized 
by only two parameters (A and (3d\/d(3). If T ~ Tj.es — 0.25^ the two level 
{F and F + 1) model is a reasonable approximation practically at any level 
distribution. In this model at np = dp the chemical potential A is between 
Ep and Ep+i. if dp = dp+i the value of A is just in the middle of this energy 
interval (A = {ep + ep+i)/2) otherwise it is closer to that level the degeneracy 
of which is lower. Decreasing np np to dp /2 makes the chemical potential and 
A + pdX/dp approach to ep i.e. in the regime dp < np < dp/2 decreasing 
np gradually shifts the minimum of ipeff{e) to ep. Subsequent decreasing 
np < df /2 leads to values of A appreciably lower than ep but \+(3d\/d(3 holds 
its position near ep. The difference between these quantities (i.e. [3d\/d[3) 
violates the left-right symmetry of ipeff{e): its left peak becomes higher than 
the right one and the more is (3d\/dj3 the more distinct is the asymmetry of 
LPeff{e). It implies that at np < dp/ 2 level F does not take part in forming 
C. If levels F and F + 1 are really far from other levels the maximum in 
C cannot practically arise as the right peak of (peff{e) disposed over £_f+i is 
small and does not give a marked amplification of C. Such a case is displayed 
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in Fig. [6l This tendency reveals itself also in Fig. [71 Interpretation of data 
in this figure can be carried out by using Fig. [8] where values of ipeff{e) are 
given for different np in Fig. [71 Fig. [7]shows also that the GCE approach (A is 
temperature dependent but pd\/dp is omitted) does not give correct results 
for all Up a.t T — > 0. Figs. [9l [10] and Figs. [TTl [12] show that for realistic 
spectra of spherical and cubic shells the forming of the local maxima involves 
not only F and F + 1 levels but also F — 1 level. 

In some cases heating mesoscopic bodies above Tres can result in a local min- 
imum in C if the preceding local maximum was created by a single electron 
level group compactly disposed near ep and far enough from other levels. Rise 
in temperature shifts therefore the peaks of (p{e) in energy regions free from 
single-electron levels that causes a decrease of C. Further heating can produce 
one more maximum in C if one or several levels appear again under the peaks 
of ip{e) as it shown in Figs. [TH [T21 These figures indicate that the second 
resonance is formed by not one pair of levels, as it takes place for the first 
resonance (5Tj.es — ep+i —£p) but in this case under the maxima of <y9(e) there 
are at least two level bunches. 



4 Shape resonances in the fermion heat capacity 

The canonical heat capacity temperature variations at T < 5 described in 
the previous section in the framework of NPM can serve as an indicator 
of the level distribution near Sp: monotonous increase of Cce is associated 
with a uniform spectrum with c? = 2, extrema of C give evidence concerning 
either local level concentration or high degeneracies of single-particle levels 
that is an inherent property of symmetric systems. The values of C are also 
affected by the particle number as the complete occupation of the Fermi level 
at T = {np = dp) in symmetric systems results in essential excess of C 
over the Sommerfeld value Cun ~ A^. Thus, due to its dependence on the 
shape of the system and nonlinearity in the particle number the mesoscopic 
heat capacity at low temperatures paradoxically differs from C of the same 
bodies at T ^ 5 where the specific heat linearly increases with T irrespective 
of the shape of the body. Examples of such shape dependence of C is given 
in Fig. [13] for 3D-systems (A^ = 40) and in Fig. [Tl for 2L)-systems (A^ = 100) 
which show the temperature variations of C for A^-particle systems of different 
but topologically equivalent shapes. 

Another example of the heat capacity shape dependence is deformation os- 
cillations in C. Before describing these oscillations from the point of view of 
NPM we will indicate that these and other low temperature oscillations in 
C are caused generally only by the quantization of the mesoscopic system 
excitation energies. 
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To develop this conception consider a mesoscopic many particle system sub- 
jected to an external impact that does not affect the volume of the system 
and thereby the particle number and that is characterized by a continuously 
varying parameter a. The latter can be an external field strength, deforma- 
tion parameter or something else. The hamiltonian or boundary conditions are 
supposed to be dependent on a in such a way that at some values of a, a = ao, 
one or several excited states of the system can merge with the ground state 
(level crossing) or, in other words, deviations of a from can give rise to a 
splitting of the ground state into two or several components (sublevels). The 
energy of the splitting Ai?sp(a) at small |a — ao| has to be much smaller than 
A£^, the energy that separates the ground state at a = ao from other highly 
excited states, Ai?sp(a) -C AE". Here the nature of the excitation spectrum 
i.e. what interactions between particles form it plays no role. The only prop- 
erty of the excited states which needs is their smooth variation determined by 
parameter a. Then, increasing a to point and further can cause at low tem- 
peratures, T < Ai?sp(a), the appearance of extrema in C . Indeed, at T ^ [\E 
and a = ao the ground state is divided by IS.E from other states. Therefore 
the heat capacity is practically equal to zero (~ exp[— Ai?/T]. Deviations of 
a from ao create a doublet or multiplet of sublevels with energies Eq ( a new 
ground state). Ex, ■ ■ ■ Ei, . . . E^, Ei — Eq < AE^p^a) = E^ — Eq {i < k) and 
the temperature transitions between these sublevels become possible. However 
two tendencies compete now: on the one hand increasing |q; — ao| ^ Ei — Eq 
intensifies the contribution of state z to C ~ {Ei — Eq)"^ , but on the other hand 
the temperature factor ~ exp[—{Ei — Eq)/T] falls with increasing {Ei — Eq). 
Consequently, at a fixed small temperature (T < AEsp{a)) this competition 
can result in the appearance of two maxima in C at values of a to the left and 
right of the point ao and the mentioned above minimum of C (in point ao) 
lying between these maxima. Thus, in a wide range of parameter a variations 
values of C will oscillate: a pair of maxima will appear every time when a 
passes through each consequent point in which the level crossing occurs. In 
Ref. PS] we showed that similar oscillations in C (and in the magnetic sus- 
ceptibility) are stimulated by the increasing magnetic field H i.e. in that case 
the role of a was performed by H. 

In NPM the level crossing {a < a^) or splitting {a > a^) in the whole 
many electron system at small deviations of a from ao are produced due to 
transformations in the single - electron spectrum: alterations in a cause first, at 
a < ao, confluence of two or several levels in a high degenerated level {a = ao) 
then, at a > ao, they split this level into levels with lesser degeneracies. 
However, the oscillations in C created by deformations are possible only near 
such points ao in which among crossing or splitted levels there are the Fermi 
level and some levels with higher energies that are free from electrons at T = 0. 
In other words, the high degenerated single-electron Fermi level arising after 
crossing at a = has to be incompletely filled at T = i.e. np < dp. 



20 



The heat capacity being plotted vs a at a fixed small temperature (T < 
AEsp{a)) varies near point like the spectral heat capacity distribution ip{e) 
vs e as on the curve of C vs a there must be, as discussed above, two maxima 
to the left and right of point in which C ~ if AEsp{a) < AE. C as 
a function of a possesses the left-right symmetry with respect to point ao 
if levels with identical degeneracies and equal spacing between them merge 
in ao and go out of this point. This symmetry implies that the amplitudes 
and distances from point ao are equal for both maxima in C. The crossing 
of two levels with different degeneracies leads to asymmetry in maxima. That 
maximum is higher which is formed by crossing the completely filled Fermi 
level {np = dp) and empty level F + 1. The higher maximum is further from 
the point ao than the second lower maximum. 

In both cases, just mentioned, the functional dependence of C on T and on 5 
(the gap between adjacent crossing levels) cannot be considered separately be- 
cause C is a function of the product [35. This quantity in the maximum {(35)m 
depends on the degeneracy of the joint level in point ao and on its occupation 
number: for crossing two doubly degenerated levels with two electrons on the 
lower of them ~ 3 as it takes place for the two level model considered 

in the previous section. Consequently the position of the maximum in C on 
the axis a depends on the temperature: if at a fixed temperature e.g. Ti its 
position is |ai — ao| then at T2 > Ti the maxima move aside from point ao 
i.e. |a2 — aol > |ai — ao| since 5 linearly varies with |a — ao| provided the 
deviations are small. Whereas the positions of the minima in C are strictly 
fastened by the positions of the points ao where (/35)o ~ /3(a — ao) = 0. This 
link of T and a implies that at each temperature (one case is e.g. given in 
Fig. [15]) the deformation |a — ao| is tuned in so that the positions of the peaks 
in the spectral heat capacity ^{s) coincide with the energies (depending on a) 
of the crossing or splitted resonance levels. Therefore if it were possible to de- 
crease the temperature in such a way that the ratio T/|a — ao| remained equal 
to a constant then maxima in C would be observed almost in the immediate 
vicinity of ao. 

To exemplify low temperature oscillations in C initiated by the deformation we 
have calculated in the framework of NPM the electron heat capacity of such 
systems as an elliptic 2Z)-oscillator system with variable ellipse semiaxes and 
3D-system in a circular cylindrical rectangular potential with variable radius 
R and height H. In both cases the spatial area or volume keeps a constant 
value at all shape variations i.e. the particle number remains invariable. 

For 2D-oscillator systems we have taken advantage of the selfconsistency con- 
dition suggested in Ref. [B]. This condition consists in that the shape of the 
particle spatial distribution (in the 2D-case it is x^a~^ -|- |/^a~^) has to be 
similar to the shape of the average potential in which particles move (for 2D- 
oscillators it is x^ujI + y'^ujy, ujq = UxOJy = const, the latter is the spatial area 
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conservation condition). Hence, the ellipse semiaxes Qx-, ay have to be pro- 
portional to reciprocal oscillator frequencies ui^, ujy, ~ cu"^, tty ~ tu"^, i.e. 
deformations of an ellipse can be defined via 2D-oscillator frequencies. There- 
fore in this case we choose the ratio of ujx/^y as a deformation parameter 
a. 

The 2Z)-oscillator eigenvalues can gain an additional degeneracy (as compared 
with the spin degeneracy) if = ^x/^y is equal to the ratio of integers z/i/z/2. 
The Fermi energy limits possible sets for ui and 1^2, besides, as discussed above, 
the heat capacity oscillations arise only near such points ao in which (at T = 0) 
the Fermi level crosses levels free from electrons. In Fig. [15] this condition is 
fulfilled only in four points of ao (6/5, 4/3, 5/3 and 2/1). In these points the 
additional degeneracy of the Fermi level is 2 for the first three values of and 
5 for cto = 2, as it shown in the bottom panel of Fig. [13 At small deviations 
of a from = 2 all 5 levels are practically equidistant. 

Near these points ao i-e- before crossing and after splitting each level is twice 
(spin) degenerated and so the heat capacity vs a has mirror symmetry with 
respect to the perpendicular to point Oq and its maxima are determined prac- 
tically only by values of (36. For the first three cases (near ao = 6/5, 4/3, 5/3) 
the values of 6, as it seen in Fig. [15] under corresponding maxima, are nearly 
identical that results in the identical heights of their maxima. In the last case 
(near ao = 2) 5 is smaller than in previous cases and the maxima are higher, 
C ~ {(36)'^ exp{—f36). The distances between maxima are explained by the re- 
lationship between S and Aa = |a — ao|. For small Aa: S ~ co'o(ao)~^^^f|Aa, 
(ao = z^i/z^2)- Since for all cases in Fig. [15] values of 6 differ insignificantly and 
the temperature is independent of a the greater distances between maxima 
correspond to the greater values of ao (and thereby the lesser values of 1^2), 
e.g. if ao = 4/3, 6 ~ 2.6u;oAa but if ao = 2/1, 6 ~ O.TwoAa. Fig. [TB] shows 
the temperature evolution of the maxima. Rise in temperature leads first to 
confluence of the low temperature deformation maxima then weak extrema 
determined by spacings between shells can be revealed and at T > 26f all 
oscillations of C vs a are smoothed. 

In the previous section in Figs. [7] and M values of C vs T are depicted for 
such symmetric systems as spheres and cubes. Now Figs. [TB] and [TH] can 
be used to explain how the variation of C vs T will change at symmetry 
breakdown. First we will consider two extreme cases: the high degenerated 
Fermi level is completely occupied at T = i.e. up = dp and the second 
case is np <^ dp. The first case is realized in that point a in Fig. [T^ where 
ao = {uJx/ujy)Q = 1 and dp = np = 12. Here small symmetry disturbance 
practically changes nothing: the maximum in C appears at the resonance 
temperature T ~ 0.2{ep^i — e^), as it observed in Fig. [16] at T/Sp = 2. In the 
second case {ux/ujy = 2, dp = 10, np = 2) in Fig. [15] at small deformations 
Aa ~ 0.05 there arise two low temperature maxima whereas at the resonance 
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temperature T ~ 0.2{eF+i — £F), where e^+i and ep are values at u^/uy = 2, a 
maximum (or even a bend) on the curve of C versus T is practically unobserved 
both in the symmetric case ujx/^^y = 2 and at small symmetry disturbance 
uix/ojy = 2 ± 0.05 (the condition np -C dp is discussed in connection with 
Fig. [6]). In the intermediate cases {np ~ dp/2) two types of local maxima 
have to occur: at T nearly equal to Tres of the symmetric system and at 
T <^ Tres corresponding to small splitting of the Fermi level at deformation. 
Thus, small symmetry disturbances do not practically change the positions 
of those maxima in C which occur in the symmetric case but at much lower 
temperatures there can appear other maxima if n^;' < dp. 

At stretching (squeezing) a circular cylinder along its symmetry axis (axis z) 
that is accompanied by decreasing (increasing) its cross-section owing to the 
volume conservation the appearance of deformation extrema in C is caused by 
accidental level degeneracies at a proper relationship between radius R and 
height H of the cylinder. Their ratio is therefore adopted to be a deformation 
parameter a = R/H. The degeneracy of single-particle levels in a potential 
with circular symmetry arises due to spin degree of freedom and in consequence 
of symmetry of the potential (conservation of angular momentum 2;-projection 
A) i.e. the degeneracy takes values equal to 4 (|A| 7^ 0) or 2 (A = 0). It means 
that levels merging in points or going out of these points can possess 
different degeneracies. As pointed out above, this difference in degeneracies 
leads to asymmetry of the maxima disposed around points ao on the heat 
capacity curve vs a. Such cases are shown in Fig. [T71 where each of the first 
three level crossing points (C ^ in these points) is bordered by two maxima 
with distinct amplitudes while the last case displays the symmetric picture. 
The following degeneracies and occupation numbers at T = can be ascribed 
to levels F and F + 1: near the first and third points ao in Fig. [T7]at crossing 
dp = 4, Up = 2, dp+i = 2, n^+i = 0, at splitting dp = np = 2, dp^i = 4, 
np+i = 0; near the second point ao at crossing dp = np = 4, dp+i = 2, 
np+i = 0, at splitting dp = np = 2, dp+i = 4, np^i = 2 and near the last 
point ao given in Fig. [T7] at crossing and splitting dp = np = dp+i = 4, 
np+i = 0. Thus, in all points (when levels F and F -|- 1 have merged) the 
high degenerated level [d = dp + dp+i) is unfilled otherwise the maxima in C 
near points ao would be impossible. The adduced numbers of d and n explain 
the values of amplitudes in maxima in accordance with the rule given above 
which can be detailed for two crossing levels F and F + 1 if one writes out the 
leading term in the equation for C: 

C ~ np{dp+i - np) [{np+i){dp - np + 1)]"^ (/5(5)^ exp(-/35), 

where 6, as before, is the difference e^+i — £f- The temperature evolution of 
the extrema in C shown in Fig. [17] develops approximately in the same manner 
as in previous case in Fig. [13 
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The examples of the heat capacity shape dependence considered in this section 
indicate that in spite of the differences in space dimensionahty and conditions 
governing the appearance of the heat capacity extrema the shape of the N 
particle system can be adjusted so (most likely it is possible to make only 
experimentally) that in a low temperature range T < ep/N the fermion heat 
capacity can attain to its maximum or any lesser value at a given (the 
needed value of C depends on the problem where these values are of particular 
importance). This property of mesoscopic systems can reveal itself also in 
granular systems. 

To show it consider two mesoscopic systems. The first is a granular system 
which is composed of close-packed metallic ultrasmall granules as it is de- 
scribed e.g. in Ref. ^5]. These granules are assumed to be in thermal contact 
i.e. they are kept at common temperature but the mutual electron exchange 
is impossible (the very weak intergranular electron coupling). Here we again 
idealize the situation and suppose that each granule contains conduction 
electrons, all they are identical in size and shape and the heat capacity of each 
granule at a fixed low temperature (< ep/N) is equal to C{N) ^ Ciin{N), 
Eq. ([1]). If the granular system consists of n such granules (the summary elec- 
tron number A^^ = nN) its heat capacity, nC{N), can be essentially distinct 
from C of the second system which is a relatively large monolithic mesoscopic 
system with the same number of conduction electrons Ng. Here the Fermi en- 
ergies of each granule and the monolithic system are assumed to be practically 
the same. Since the mean level spacing in the monolithic system is n times 
smaller than in the granule the fixed temperature can appear so high for the 
monolithic system at n ^ 1 that its heat capacity can take only the shape- 
independent Sommerfeld's value Cun^Ng) = nCuni^N). Therefore the ratio of 
the electron heat capacities of the granular system and the monolithic one is 
equal to the quantity C{N) /Cun^N) given in some figures of this and previous 
sections such as Figs. O [9l [TT|, [121 Thus these figures gain the second content: 
they give not only values of C/Cun vs T for a system with a fixed particle 
number but the same curves describe the ratios of C in granular and mono- 
lithic systems, where the former are composed of small A^ particle granules: 
briquettes (Figs. [5| [T2l). cubes (Fig. [9]) and spheres (Fig. [11]). Additionally in 
Fig. [18] we show the variety of values C /Cun for cubic systems in a wide range 
of AT < 3 • lOl 



5 Conclusion 

For description of fermion properties of mesoscopic A^-particle systems which 
are in thermal equilibrium with surroundings (the canonical ensemble) in the 
framework of NPM we have suggested a new formalism representing the 
single-particle level occupations as a continuous function of energy e similar 
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to the Fcrmi-Dirac distribution. In NPM many static thermodynamic quan- 
tities can be obtained by applying this representation. In particular it allows 
the electron heat capacity C to be calculated as a convolution of the state 
density with function 97(e), the spectral distribution of C that gives the con- 
tribution of each single-electron state to C. The two hamped shape of ip{e) vs 
e gives possibility to account for the appearance of local maxima in C within 
the low temperature range (T < dp-, the latter is the mean level spacing near 
Ep) by virtue of the resonance amplification of contributions of states nearest 
£i? if the positions of two peaks in ip{e) coincide with energies of either two 
high degenerated levels near Sp (or level bunchings) or even two solitary dou- 
ble degenerated levels. However C vs T for the uniform spectrum with equal 
spacing between spin degenerated levels does not reveal noticeable maxima. 
It is worth mentioning that local maxima in C are observed more distinctly 
on curves C/Cun (or C/T) vs T where C^n is the Sommerfeld's electron heat 
capacity linear in T. 

The temperature of the maximum in C (the resonance temperature T^es) 
makes up a part of Ae, the difference between energy levels in the vicinity of Sp 
and depends on ni;', the Fermi level occupation number at T = 0. So if = cii? 
is high {dp > 2) {dp is the Fermi level degeneracy) then Tres — ^■'^{^f+i -~£f) 
whereas for the case of two solitary spin degenerated levels {cp and £p+i) 
with Up — 1 this temperature is equal to ~ 0.4(£i;'+i — Sp). Thus, establishing 
Tres and studying A'"-dependence of C in samples of identical shapes can give 
information concerning the single-electron level structure near £p. 

Within the low temperature range the heat capacity appears to be shape- 
dependent at an invariable spatial area (for 2Z}-systcms ) or volume (for 3D- 
systems). At continuous alteration of some parameter a characterizing defor- 
mation of the system, C can display oscillations vs a with minimum values of 
C at such deformations where the system gains an additional symmetry. 

It is shown that at T < Sp/N the ratio of the electron heat capacities of 
a granular system consisting of n identical insulted iV-electron subsystems 
{n ^ 1) and a monolithic system with the same summary electron number 
{n ■ N) is equal to C{N)/Ciin{N) which can be stretched from zero to several 
units. 

On the whole, the presented material can be viewed as evidence of possi- 
ble essential deviations of low temperature (T < 8p) properties of electron 

mesoscopic systems from those established for T > Sp that is caused by the 
irregularities in single-electron spectra and this has to be taken into account 
at consideration of both the heat capacity and all other thermodynamic prop- 
erties of finite electron systems. 
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Fig. 1. The effective grand canonical and canonical heat capacity (a) and CjCun 
(b) of N-even system vs T in a model with the equal level spacing 5, the level 
degeneracy of each level d = 2. np is the occupation number of the Fermi level at 
r = 0. 
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Fig. 2. The spectral heat capacity in the EGCE and CE at different temperatures 
for the equal level spacing model. Vertical lines mark the level positions. 
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Fig. 3. Top panel: The heat capacity in the EGCE and CE vs T in tlie two level 
model with np = '2, and d = 2. Vertical lines mark the resonance temperatures. 
Bottom panel: The spectral heat capacity in the EGCE and CE for resonance 
temperatures. 
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Fig. 4. Top panel: The heat capacity in the EGCE and CE vs T in the two level 
model with np = 1 and d = 2. Vertical lines mark the resonance temperatures. 
Bottom panel: The spectral heat capacity in the EGCE and CE for resonance 
temperatures. 
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Fig. 5. Top panel: The heat capacity (in units of Cun) of A^-even and A"-odd systems 
vs T/dp for a briquette. The lateral lengths {L^ : Ly : = 1 : 0.6e : tt) are chosen so 
to provide only the spin level degeneration, Sp = 4£i?/3iV. Bottom panel: Fragment 
of the single-particle spectrum of the N = 1440 briquette. 
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Fig. 6. Impact of the Fermi level filling on the low temperature resonance in the 
canonical heat capacity. Two level model, dp = dp+i = 10. 
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Fig. 7. The same as in Fig.[6]but for the = 40 shell in the spherical cavity, n^? is the 
particle number at T = on the Fermi level with degeneracy dp = 6. For a sphere 
with = 40 the mean level spacing Spisph) = sp/N^^'^ ~ £f+i — £f- Results of the 
canonical {CE), equivalent grand canonical {EGCE) and grand canonical {GCE in 
which A is temperature dependent but f5d\/df5 is omitted) methods are displayed. 
It is obvious that the GCE method gives correct values of C at T — > only at 
np = dp and np = dp/2. In both other methods C — > at T — > 0. 
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Fig. 9. Impact of the Fermi level filling on the canonieal heat capacity (in units 
of Clin) for the sheU of a cube with = 991 ^ 1008. Spicuhe) = 20eF/3iV is 
adopted as the mean level spacing near the Fermi shell. For a cube with N = 1000 
5F{cuhe) ep+i — sp- Cegce and Cce practically coincide since ^ 1. 
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Fig. 10. The spectral distribution of the effective grand canonical heat capacity 
iPf.ff{e) calculated for = 1;9 and 18 (see Fig. [9]). 
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N = 1 32 




Fig. 11. Top panel: The heat capacity (in units of Cun) with two resonances for 132 
fermions in spherical cavity, ti and t2 are the positions of maxima. Bottom panel: 
The spectral heat capacity ipeff(e,t) at these temperatures. The first maximum 
is practically at ti = {ep+i — eF)/5SF{sph) but t2 corresponds to the presence of 
two groups of levels in the vicinity of the maxima in ip^ff{£,t). Vertical lines mark 
the single-particle levels, their heights are proportional to the level degeneracies. 
Spisph) = Ep/N^^'^ is the mean level spacing between spherical shells. 
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. 12. The same as in Fig. [U] but for the N = 10^ briquette, dp = Aep/^N 
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Fig. 13. Top panel: The heat capacities of = 40 systems in 3-D-cavities of differ- 
ent shapes (briquette, cube, cylinder and sphere ). The dashed vertical hues mark 
values of t = T/5f = 0.2(ei?+i — eF)/6p, Eq. (fSTj) . These temperatures practi- 
cally coincide with the positions of the corresponding maxima in C. In this fig- 
ure dp = Aep/'iN. Bottom panel: Fragments of the single-particle level schemes 
of sphere, cube, cylinder (the diameter is equaled to the height) and briquette 
{L^:Ly:L, = l: 1.1 : 1.2). 
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Fig. 14. The heat capacities of A'^ = 100 systems in 2D-cavitics of different shapes 
(circle, square, rectangle with Lx : Ly = 1 : 0.318). dp = 2eF/N. 
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Fig. 15. (a) - Shape oscillations of the low temperature canonical heat capacity of 
the system with iV = 42 fermions moving in 2D oscillator potential, ujx^ are 
the oscillator frequencies, dp = ep/^N. (b) - Electron level energies near the Fermi 
level (bold line) vs uj^l^y Solid circles are the crossing points of the Fermi level 
and upper levels free from electrons at T = 0. 
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Fig. 16. The temperature evolution of the shape oscillations of C for the same 
system as in Fig. [15l ujx and ujy are the oscillator frequencies. 
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Fig. 17. Shape oscillations of the low temperature canonical heat capacity vs R/H 
for the N = system in the 3D rectangular circular symmetric cylindrical cavity 
with radius R and height H. Sp = 4ei?/3iV. 
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Fig. 18. N- oscillations of the canonical heat capacity {C/Cun) in cubes, N < 3-10 . 
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